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Abstract. We study an indefinite Sturm-Liouville problem due to Richardson 
whose complicated eigenvalue dependence on a parameter has been a puzzle for 
decades. In atomic physics a process exists that inverts the usual Schrodingcr 
situation of an energy eigenvalue depending on a coupling parameter into the so- 
called Sturmian problem where the coupling parameter becomes the eigenvalue 
which then depends on the energy. We observe that the Richardson equation is of 
the Sturmian type. This means that the Richardson and its related Schrodingcr 
eigenvalue functions are inverses of each other and that the Richardson spectrum 
is therefore no longer a puzzle. 
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1. Introduction 

In the early years of the twentieth century an eigenvalue problem based on an indefinite 
second order differential equation of the Sturm-Liouville(SL) type was introduced by 
the mathematician R. G. D. Richardson [T]. The history of this early period has 
been reviewed by Mingarelli [2] with later developments summarized by Binding and 
Volkmer [3] . The spectrum of the Richardson problem depends on a real parameter in 
a very complicated way that has never been fully understood and has even been called 
amazing [3] . Our goal here is to provide a non-rigorous explanation of this spectrum. 

Although the Richardson problem predates quantum mechanics, we make a 
quantum connection by asserting that, in the language of atomic physics, it is a 
Sturmian problem 5J. In such a view the Richardson eigenvalues are coupling 
parameters of a potential and the parameter that causes the complication is an energy. 
In the atomic physics of the Coulomb potential one is motivated to go from Schrodinger 
to Sturmian for computational advantage. Here we reverse the direction and go from 
Richardson as a Sturmian problem to Richardson as a Schrodinger problem where 
now more naturally the parameter is the coupling and the energy is the eigenvalue. It 
is then seen that the Richardson and Schrodinger eigenvalue functions are inverses of 
each other and thus both contain the same information. We are then led to understand 
Richardson in terms of Schrodinger. 

In the Richardson problem the independent variable is the energy E and we call 
the dependent ones the eigencouplings X n (E). In this case for E real and fixed we 
have one energy and many potentials while in the Schrodinger problem with A fixed 
we have one potential and many energies E n (X). We use the term spectral inversion 
to mean going from one description to the other. The complications of the Richardson 
spectrum [4j involve short segments in E separated by square root branch points in 
which X n (E) is alternatively real, imaginary or complex, and these possibilities recur 
over and over as E varies through real values. In the related Schrodinger problem, 
E n (X) for real A is much simpler and this difference in complication for real values of 
the parameters must stem from a difference in the two SL equations that generate the 
spectra. In SL terminology [6] the Schrodinger equation is intrinsically right-definite 
and if certain technical requirements are satisfied, the E n (X) are real for real coupling. 
In the Richardson case the eigencoupling multiplies the spacial part of the potential 
and if this function changes sign as it does in the Richardson problem, the SL equation 
can be both left- and right-indefinite and this leads to a spectrum that is both non-real 
and singular as E varies through real values. In this paper we argue that one of the 
two descriptions may not be simpler than the other because inverse functions contain 
the same information. This equality implies that the Schrodinger version must be 
suitably complicated in the complex domain of A where the E n (X) are less tightly 
constrained. In this way a balance of complication can be maintained. It must be 
so that the Schrodinger complication is more subtle and is not as evident in its SL 
equation. 

To examine these questions we need to study the eigenquantities in the complex 
domains of their independent variables, or more generally on the Riemman surfaces 
on which they take their values. The analytic properties are important and we need 
to focus on their square root branch points and critical points since one type of point 
becomes the other under spectral inversion. It will turn out that critical points are the 
sources of the Schrodinger complication while square root branch points are responsible 
for the complication in the Richardson problem. Thus spectral inversion relates one 
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source of complication to the other. From the mathematical point of view it also 
relates an indefinite SL problem to a definite one. 

In section [5] we relate the critical and square root branch point of the Richardson 
and Schrodinger problems in terms of single and double zeros of a certain spectral 
function. It is also pointed out that Schrodinger critical points for real A lead to 
curves in the upper and lower planes of A on which the eigenvalues are real. Section [3] 
discusses the Richardson and Schrodinger differential equations. In section |4] we treat 
the problems more generally and develop equations for the generalized eigenfunctions 
that are present at branch points. The Schrodinger spectrum is studied in detail in 
section [5] where branch points and critical points are located and the Riemann surface 
structure is discussed. Section [6] treats the analogous properties of the Richardson 
spectrum and the structure of its Riemann surface. A qualitative demonstration that 
the inversion of the Richardson A„ (E) leads to the Schrodinger E n (A) is presented in 
section [7] Some questions concerning oscillation counts and defective multiplicities 
are treated in section [5) Things are summarized in section [5] 

2. Level crossing and its analog 

Our goal in the section is to show the relation between singular points and critical 
points in the Schrodinger and Sturmian descriptions. In the mathematical literature 
there is some discussion of these questions at an abstract [7] as well as at a more 
concrete [3] level. 

We start by considering a function A (A, E) of the coupling parameter A and 
energy E whose vanishing indicates a discrete level of the system. For example, this 
function could be the secular determinant formed in an attempted diagonalization of 
the Hamiltonian in some basis. The vanishing allows us to define an energy E{\) if 
A is independent or the eigencoupling X(E) if E is independent. To be specific in the 
Schrodinger case suppose that values Ao and Eq cause A to vanish. Expanding near 
the zero we have 



with the subscripts indicating partial derivatives. Choosing A so that A(A, E) vanishes 
allows a local solution for -E'(A) 



Level crossing at Ao means that A has a double root so that Ae also vanishes which 
makes -E'(A) singular. This argument does not disclose the nature of the singularity 
but generically they are square roots which we will assume to be the case since we 
have no numerical evidence to the contrary. 

If the above steps are repeated for the Sturmian problem and if A(Ai,-Ei) = 0, 
we obtain 



A(A,£) = A A (A-A ) + A £ (i; 



E ) + ■ ■ ■ , 



(1) 



E(X) =£(A )--A(A-Ao) + - 



£(Ao) + -tv (A-A ) + ---.(2) 
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Coupling crossing occurs if A\ = which makes X(E) singular. The possibility of 
coincident zeros of 1% and I\ has been discussed and eliminated at least in one case 



We now return to the Schrodinger problem and ask what property of E(X) at 
a regular point of A implies X(E) is singular at the image point. Returning to ([2]) 
we invoke Ae / with A a = which gives E'(X) = 0. Thus a critical point of 
E(X) implies that X(E) is singular. Of course a critical point position Ao and critical 
value -Eo(Ao) in the Schrodinger description also gives the branch point position Eo 
and singular value Xq(Eq) in the other description. In a similar way in ([3]), a critical 
point of X(E) signals level crossing in E(X). We assume the absence of triple or higher 
crossing in a problem with one coupling parameter present which implies that the 
critical points are non-degenerate so that second derivatives are non-vanishing at the 
points. In what follows we will sec that the cigenquanties in the two problems take their 
values on Riemann surfaces so that when we go from one to the other a (regular non- 
critical, regular critical, singular) point will become a (regular non-critical, singular, 
regular critical) point on the other. We will refer to the Schrodinger Riemann surface 
as R$ an d to the coupling one of X n (E) as Rq- 

We have stated that critical points are the sources of Schrodinger complication 
and this arises from the leading quadratic behavior in their vicinity. This induces 
a more rapid angular dependence and causes contours of Re E n (X) = constant and 
of Im E n (X) = constant to intersect in a double-x fashion. If E n (X) is real and has 
a critical point for A real, it acts as a source of secondary curves into the upper 
and lower planes on which E n (X) is also real. The real contour and the real axis are 
locally orthogonal but the curve then bends away depending on local details. Problems 
exhibiting PT-symmetry [9] usually have potentials with an even-odd spacial structure 
that leads to E n (— A) = E n (X). This makes A = a critical point but also forces all 
higher odd derivatives to vanish causing curves into the upper and lower plane to 
coincide with the imaginary axis over an extended interval. It is this special type of 
critical point that gives rise to the properties of PT-symmetry. If any odd derivative 
higher than the first is non- vanishing, then the curves will not coincide with the axis 
and the eigenvalues will be non-real on it. If one is interested in real E n (X) in the 
complex domain of A where the Schrodinger operator is not self-adjoint, critical points 
of -En (A) for real A are of special importance. Under spectral inversion these points 
will become square root singularities for E real and thus contribute to making X n {E) 
complicated. 

3. Richardson and Schrodinger problems 

The Richardson problem is classed as a right-indefinite SL equation and may be written 



0. 



(t>'n(E,x) -E(f> n (E,x) = X n (E) sgn(x)(j) n (E,x), 



(4) 





£ is a parameter that is usually real and the function sgn(a;) changes sign once in 
the interval and is termed an indefinite weight. Its presence leads to the complication 
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of the Richardson spectrum; if it is replaced by 1, the problem is elementary. The 
eigencouplings X n (E) are the quantities of primary interest. 

When we switch to the Schrodinger description, the SL convention is to put the 
eigenvalue term on the right so that E — > E n (X),X n (E) — > A, and we obtain 



The indefinite weight has become incorporated into the potential V(x) = — X sgn(a;) 
and may seem to be of less significance; no matter how many times V(x) changes sign, 
([5]) is right-definite and is much better behaved than ((4|). This potential for purely 
imaginary coupling has been independently suggested as a problem exhibiting PT- 
symmetry by Znojil [10j . followed by joint work with Levai [11] and we can compare 
our results to theirs in Section [5] 

We now make a number of remarks concerning spectral inversion. If we go 
from Richardson to Schrodinger we have X n (E) — > E n (X) so that place — » value and 
value — > place. The Richardson X n (E) represent complicated values and real places 
which when inverted become real values and complicated places. This Schrodinger 
combination has as its source the critical points for real A that produce curves in the 
upper and lower plane on which E n (X) is real. It follows that we view the Richardson 
spectrum, we are viewing the complex geography of these Schrodinger curves. 

Suppose we have determined that particular values of Ao and Eq are image 
points on the two Riemann surfaces. If we insert these numbers into equations (j4|) 
and ([5]), the left-right details are irrelevant and the equations coincide so that the 
eigenfunctions can be taken to be identical. In the real domains of A and E if we let 
the oscillation count denote the number of internal nodes of either eigenfunction, then 
that count as determined at image points is invariant under spectral inversion, as is 
the eigenfunction. 

4. Differential equation 

We need to examine the differential equations of the two problems in more detail so 
that the role of the critical and branch points can be made clear. We also study the 
generalized eigenfunctions present at branch points where either spectrum is defective. 
In the Schrodinger case Simon [8] has shown in a study of the anharmonic oscillator 
that the inhomogeneous equation for the generalized eigenfunction can be obtained 
by taking the derivative of the Schrodinger equation with respect to E. In both 
Schrodinger and Richardson cases we will use procedures analogous to this. 

A more general potential V(x) — Xv(x) + w(x) is chosen where the two functions 
are real and A may be complex. Initially, we leave open whether A or E is independent. 
In abbreviated notation the Schrodinger equation is 



- ip n (X,x) - X sgn(x)ip n (X, x) = E n (X)ip n {X,x), 
-1<x<1 ^„(A,±1) = . 



(5) 



- ip (x) + Xv(x)ip(x) + w(x)ip(x) — Eip(x) = 0, 

= ijj{b) = 0, — co <a<x<b<oo . 



(6) 
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In addition to this it is useful to have the following result. Suppose a second function 
£(x) obeys an equation that is identical to © on the left but has an inhomogeneous 
term F(x) 

- {"(x) + Xo(x)£(x) + w{x)£(x) - E£(x) = F(x), (7) 
£(a) = £(&) = 0. 

We want a solution to both equations for the same values of A and E. If we multiply 
(JHJ) by £(x) and ([7]) by ip(x), subtract and integrate over x from a to b we obtain 

(F\ip) R = I dxF(x)ip(x) = 0. (8) 

J a 

In other words for a pair of solutions to exist, the inhomogeneous term must be 
orthogonal to the homogeneous solution with respect to a real (R) inner product. 

We now go to a point E on M-c and wish to study the derivative 8e^{x) = iPe(x). 
Letting D = —d% + Xv(x) + w(x) — E and acting on © with 8e we get the pair 

DV(x) =0, (9) 
T)^e(x) = -(dX/dE)v(x)ij(x)+ip(x). 

At a critical point the derivative on the right vanishes giving the simpler pair 



Dip(x) = 0, (10) 
D^e(x) = ip(x). 



Solution of these equations would produce four things: the critical point position E 
on M-c, the critical value X(E), as well as ip(x) and iI>e{x). The image branch point 
on Kg has the same values of A and E, dX/dE vanishes there as well, and (fTU)) has 
the standard form [8] for a eigenfunction-generalized eigenfunction pair at a branch 
point in a Schrodinger problem. We then associate this iPe(x) with the generalized 
eigenfunction on R<?. Putting in a full notation and renaming %j)E — * ip, the pair (|10|) 
becomes 



-ip'n(X,x) + Xv(x)ip n (X,x) +w(x)t/j n (X,x) = E n (X)4> n (X,x), 

-ipn (A, x) + Xv(x)ip n (X, x) + w(x)ip n (X, x) - E n (X)if) n (X, x) = tp n (X, x). 
The constraint of |5]) gives 

(V»n(A)|^„(A)) fl = / dxip n (\x)i) n {X,x) = 0. (12) 
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We can carry out an analogous procedure starting at a critical point A on R5. We 
write ip\( x ) = 9\ip(x) and apply the derivative to ^ which gives 

£>V>0) =0, (13) 
S^a(x) = {dE/dX)ip{x) - v{x)ij){x). 



Whether at the critical point on R5 or the image branch point on Rp, dE/dX = 
and the simpler pair is 

S)ip(x) = 0, (14) 
T>ijj\(x) = —v(x)ip(x). 



As before, ip\(x) has the interpretation as the generalized eigenfunction on Rc. We 
rename if) — > <j> and ip\ — > to give the pair of equations 



(15) 



-4>' n (E,x) + w(x)<jt n (E,x) - E<fi n (E,x) = -\ n {E)v{x)4> n {E,x), 

-4> n {E,x) + w(x)4> n (E, x) - E<j) n {E,x) + X n (E)v(x)4> n (E, x) 

= -v(x)cj> n {E,x). 

Imposing the constraint of © gives 

{(j> n {E)\v\(j) n {E)) R = [ dx<p n {E,x)v(x)<j> n (E,x)=0. (16) 



The quantity ip\(x) at a point A is the first order correction to ip(x) in Rayleigh- 
Schrodinger perturbation theory. If one calculates this at a critical point on Rg by 
the usual procedures, one has wittingly or unwittingly also calculated the generalized 
eigenfunction at the image branch point on Rc- 

For non-critical and non-singular points in cither problem, the second equation 
in each pair is ignored and we deal only with the first of (fTTj) or the first of (fl5|) . 
The Schrodinger version is intrinsically right-definite while the Sturmian one is right- 
indefinite if v(x) changes sign. 

Although we cannot write expressions for E n (X) or X n (E), we can obtain useful 
ones for their derivatives by the following procedure. We return to the pair ([13"]) valid 
at a point A on R5. The inhomogeneous term must conform to ([8]) which permits a 
solution for the derivative 

dE n (X) _ ftMA)|z# n (A)) fl 

dX ftMA)|Vn(A)) fl ' 1 ' 

Going to the other pair (|10|) and performing a similar step gives 
dX n (E) (K(E)\ME))b. 



dE (ME)\v\ME))i 



(18) 
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If A and E are image points, we can take |^„(A)) = \<j) n {E)) and the two expressions 
say simply that E'{\) = \/\'(E). The integrals implicit in these expressions may 
not themselves be singular since the integrands are entire and the intervals are 
finite. Singularities must come from vanishing denominators which arise from the 
orthogonality conditions in (|T2]) and (fT6|) . In a similar way critical points arise from 
vanishing numerators. Although we do not use these derivative expressions in practice, 
they contain the essence of the relationship between the Schrodinger and Richardson 
problems. 

If A is real, the first of (fTTj) is right-definite and E n (X) and \ip n (X)} are real also 
so that the denominator in (| 1 T[) will not vanish and E n (X) will be non-singular as 
expected. According to (fT2|) the denominator will vanish at branch points which must 
then occur only off the real axis. If A is real and v{x) is of one sign, the numerator 
cannot vanish which excludes critical points and leaves E n (X) both non-singular and 
non-oscillatory. If A is real and v(x) changes sign, the numerator may vanish leading 
to curves in the upper and lower plane on which E n (X) is real. All such points on R5, 
whether A be real, imaginary or complex, must map to the real axis of Rc on some 
sheet and thereby generate X n (E) that match the three possibilities of A on M5. This 
is how the Schrodinger complication leads to the Richardson one. 

If we attempt a similar analysis of (fl8|) for real E, we can draw no general 
conclusion because the first equation of (fT5*)l may be sufficiently indefinite to give 
complex results. If v(x) is one sign, right-definiteness should make everything real with 
numerator and denominator non- vanishing in (|18|) . This gives again a non-singular 
and non-oscillatory function and it is natural that two such featureless functions can 
be inverses of each other. Suppose for E real that v(x) changes sign so that the 
denominator vanishes giving a branch point on Wc- The vanishing denominator here 
at the point E is the same as the vanishing numerator in (fT7|) at the critical point 
A. This is just another way of saying that the eigenquantities are the inversions of 
each other. We conclude that sign changes in v{x) can induce both Schrodinger and 
Richardson complications. 

There is another tool that is useful in characterizing the spectrum in each 
problem. If v(x) is of one sign, then each eigenquantity is a Herglotz fucntion [5] 
of its independent variable. To study this we start with the first equation of (fTTj) and 
act on it with Im f dxip^X, x). The integral of the derivative term and of the w{x) 
term are real and do not contribute, leading to 

\mE n (X) = \mX — , (19) 

{lp n {X)\lp n {X))H 

where the H denotes Hilbert space inner product. If v(x) is positive in (a, b), the 
ratio is real and positive for all n and all complex A. This makes E n (X) a Herglotz 
function and it can be real only on its real axis. Critical points would produce curves 
that violate this and are not allowed. The Herglotz property and PT-symmetry are 
therefore incompatible. If v{x) was always negative, we could come to the same 
conclusions and refer to the function as anti-Herglotz. 

To study A„ (E) we start from the first equation of (|15p and take similar steps to 

get 



(20) 
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If v(x) is always positive then X n {E) is also Herglotz and so this property is invariant 
under spectral inversion. In other words the three possibilities of upper plane, real 
axis, and lower plane map, respectively, into the same possibilities. 

The right-definite arguments that we have used lead to real eigenquantites when 
their parameter is real. If either of the problems had real curves off the real axis of 
their parameter, the complex places and real values would invert to real places and 
complex values. Since this is impossible, both arguments rule out such curves. This 
argument can also be applied to the Richardson problem to exclude the existence of 
curves off the real axis of E on which lmX n (E) = 0. The spectral inversion of this 
would be unacceptable to the Schrodinger problem. 

Suppose v(x) does change sign so that critical points generate real eigenvalues 
for non-real A. E n (X) is then not Herglotz but (|T9|) is still valid and it follows that 
(ipn(X)\v\i(} n (X)) h = at every point on such curves. The algebraic sign of this 
quantity is called the signature [Hj and can be used to obtain oscillation results for 
indefinite problems. 

An example of a problem that is both right- and left-definite in its Schrodinger 
and Sturmian versions is a quartic anharmonic oscillator [13j with V(x) = Xx 2 + x 4 . 
In this problem on x G [— oo, oo], the additional function w(x) = x serves to stabilize 
the spectrum at A = 0. Each eigenquantity is a Herglotz function of its parameter 
and is therefore non-oscillatory and non-singular in the real domains. 

5. Schrodinger spectrum 

In this section we focus on the eigenvalues of the Schrodinger equation. Critical and 
branch points are located and their properties studied. The parameter E that is real in 
Richardson becomes the real E n (X) in Schrodinger. We must then search in complex 
A for real E n (X). The Riemann surface of E n (X) is also described. 

Given the equation InxE n (A) = 0, we ask for what complex A is this true? Since 
ImE^A) is itself real, we make the replacements lmE n — > F, ReA — > x, ImA — > y, 
giving the real equation F(x,y) = 0. The implicit function theorem implies that the 
solution is a curve and in the complex case we call such a curve a real locus and allow 
it to consist of the following possibilities: the real axis A, the imaginary axis B, or 
an actual curve C not on either axis. This is not to say that all problems have the 
three types but only that these are the possibilities. A self-adjoint problem with real 
coupling will have A's while a problem exhibiting PT-symmetry may have B's. We 
know of no discussion of possibility C but this type must be present in the Schrodinger 
problem if the Richardson spectrum is to be understood. Since it will be argued that 
the Riemann surface of E n {X) has an infinite number of sheets, we should expect and 
infinite number of each of the three types. 

It is natural that an A and a B will intersect at A = and our discussion in section 
[2] implies that all odd derivatives of E n (X) will vanish there if E n {— A) = E n {X). If 
some odd derivative other than the first is non- vanishing, then a C would arise at 
A = and not a B. In the Schrodinger problem, the C's have as their sources the 
critical points on the real axis for A ^ 0. 

The Richardson spectrum starts at E — with E increasing through real values. 
The Schrodinger image of this must have E n (X) real and therefore involve movement 
on the A's, B's and C's in such a way that E n (X) always increases. This Schrodinger 
movement in place (A) is the encoding of the Richardson complication in value (X n (E)). 
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In sections [3] and 0] we have discussed the relevant differential equations in detail 
but for the Schrodinger problem the potential V{x) = —A sgn(x) is piecewise constant 
so that the wave function can be written down explicitly. This is the approach of 
Atkinson and Jabon [1] which we adopt. In the Schrodinger problem for fixed A we 
need a function A(A, E) that vanishes at an eigenvalue and this function is supplied 
by the difference of the logarithmic derivatives of the left and right eigenfunctions 
at x = Q. For A fixed we search for zeros in E and for E fixed we look for zeros 
in A. These are simple zeros but for branch points in either problem the algebraic 
multiplicity is 2 and we need double roots. In Schrodinger this requires A = Ag = 
while in Richardson we need A = A a = 0. There are numerical ways [13j of searching 
for such points while having no prior estimate of their location. The above method has 
problems near points where the logarithmic derivatives are singular. And alternative 
scheme is to expand the wave function in solutions of the problem when A = 0. The 
new A(A, E) is then the determinant formed from the characteristic equation. 

Given the potential V(x) = —A sgn(x), we have solved numerically for E n (X) for 
n = 1 — 4 with results plotted in figure [T] As |A| — > oo, one side or the other of the 
well becomes infinitely deep so that E n (X) — > — oo for all n. These E n (X) are referred 
to as eigencurves by Binding and Volkmer [3] who have shown that the n-th curve has 
n local maxima. We see numerically that there are n — 1 local minima and so it is 
assumed that E n (X) has 2n — 1 critical points. In [5] it is shown that the position of 
the maxima in the (A, E) plane form a rectangular lattice so that these positions and 
values are known exactly. This means under spectral inversion that the information 
for the Richardson branch points is also known. The minima do not form such a 
lattice and their positions and values have been determined numerically. 
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Figure 1. E n (X) versus A for n = 1 — 4. 
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Binding and Volkmer [3] have argued and to some extent proved that there are 
complex connections between the eigencurves of figure [1] We wish to elaborate on 
this and argue that all of the eigencurves are different branches of the same function 
and that the complex connections are facilitated by square root branch points off the 
real axis of A. In order to form a Riemann surface for E n (\) we make many copies 
of the A plane, label each sheet by n, place the n-th eigencurve over the real axis of 
the n-th sheet, and then search for square root branch points that couple one sheet 
to another. After a lengthy numerical calculation we find the branch points shown 
in figure [5] for |A| and n not too large. The two integers indicate the sheets that are 
joined at the square root. E n (X) is subject to two constraints: since E n (\) is real on 
the real axis we have E n (X*) = E^(\), and in [3j it is proved that E n {— A) = E n (X). 
These conditions require that branch points on the imaginary axis appear as mirror 
pairs and also that those in the first quadrant at A = x + iy must appear at the four 
positions A = ±x ± iy. With everything in place there are An — 2 branch points on 
sheet n. The singularities tend to move out with n and become more numerous so it 
is possible that they accumulate at infinity on the asymptotic sheet. 
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Figure 2. Positions of branch points of E n {\) for ReA > and ImA > 0. 

On relevant sheets we imagine branch cuts starting at each square root and 
extending radially to infinity. The process of level crossing involves starting on the 
real axis of sheet n with oscillation count n — 1, encircling a single branch point of 
the type n — n + 1, changing sheets, and then returning to the real axis of sheet n + 1 
with oscillation count n. This process will play an important role in understanding 
the Richardson spectrum. 

Some numerical values for the low-lying branch point positions and values are 
listed in table [1] for those on the positive imaginary axis and in the first quadrant. 
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Table 1. Low-lying branch point for ReA > and ImA > 0. 



n-n+1 


RcA 


ImA 


KcE 


ImE 


1-2 


0.0 


4.475309 


6.401903 


0.0 


2-3 


9.264139 


6.834853 


17.617719 


0.960866 


3-4 


0.0 


12.801544 


30.979714 


0.0 


3-4 


28.204239 


7.318111 


37.550337 


1.481781 


4-5 


16.798312 


15.527134 


52.144783 


1.436416 


4-5 


57.481587 


7.358543 


67.167957 


1.676906 



Using different methods, Znojil and Levai [10l QT] have also determined some of these 
values and our results are in good agreement. 

That E n {— A) = E n (X) implies that all odd derivatives vanish at A = on every 
sheet. This means that a real locus of type B is present on the positive and negative 
imaginary axes until mirror branch points are reached. We can read off from the figure 
when that happens for sheets n — n+1. E n (A) is complex above and below the branch 
points so that the real locus includes but then stops and these points. 

The presence of critical points on the real axis for A ^ implies that real loci 
of the type C move into the upper and lower planes. Such points are absent on the 
n = 1 sheet and we find that the next two sheets, n — 2 and 3, contain a pair of real 
loci that are closed curves as shown in figure [3J The symmetries of E n (A) imply that 
the curves appear as left-right mirror images. The four n = 2 — 3 branch points with 
their radial cuts facilitate closure. These curves constitute fundamental cycles on the 
Riemann surface on which the eigenvalues are real. 

Suppose we are on A or a B so that E n (X) is real. The property E n {— A) = E n (X) 
implies that there are two values of A with the same real energy. In Richardson this 
means that there are two real or two imaginary eigencouplings on different sheets for 
some real E. On a real locus of type C things are different. The end points of the two 
lines through the origin in figure [3] represent positions on the two curves with couplings 
±A and ±A*. The symmetries of E n (X) imply that all four couplings have the same 
real eigenvalue. In Richardson we will then have four distinct complex couplings on 
different sheets for some real E. This argument works for any point on the C with 
the exception of the real axis. We define these to be ^-points so that only doubling 
occurs. To summarize, the A's, -B's and C's have reflection doubling but for the C"s 
it becomes quadrupling because of ±A*. 

Since we have assumed that the number of critical points on sheet n is 2n — 1, it 
is expected that the number of real loci will increase on the higher sheets. In figure 2] 
are shown the positions of all critical points for A > and n = 1 — 6. The pairs that 
form a real locus are shown joined by an arrow. We do not display these curves but 
they are similar to the ones shown in figure [3] and branch points from the previous 
discussion produce the necessary level crossings that cause the curves to be closed. 

There is an important difference between the locus types which is worth pointing 
out. If we go to the n = 1 sheet for some real Ao and look up to all sheets for that 
coupling, an infinite number of real eigenvalues will be present. If we repeat on the 
imaginary axis and if |Ao| is not too large, the number will also be infinite. This is 
somewhat indicated in figure [4] where the critical points at A = are present vertically 
on all sheets. If we ask for the number for the C's , a look at figure 2] shows that 
the critical points are not vertical in a similar way and a given C does not appear on 
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Figure 3. Real loci of type C on sheets n - 
points from which cuts go to infinity. 



2 and 3. The open circles are branch 



other sheets. This implies that no point on a C has an infinite number that is similar 
to the other cases. To summarize, there are an infinite number of real eigenvalues for 
all j4-points, some i?-points, and no C-points. 

We have argued that there are 2n — 1 critical points of E n (A) for A real. If there 
were additional cases off the real axis, these would invert to become branch points 
of X n (E) for complex E which could lead to coupling crossing. Hundreds of searches 
using the methods described in [13] have been carried out which yield only critical 
points for real A. This implies that Richardson branch points occur only for real 
values of E. Level crossing branch points for A and E n (X) complex will invert to 
critical points in the complex domain of E. Since they play no important role in the 
Richardson problem, we have not studied them. 



6. Richardson spectrum 

Here we describe the behavior of the Richardson eigencouplings as functions of the 
parameter E and suggest a suitable Riemann surface for the problem. Some numerical 
results are also presented for the eigencoupling on a few sheets. 

Starting with the Schrodinger eigencurves of figure [1] we proceed to rotate the 
graph and adjust the labels to get the Richardson eigencurves of figure In 
Schrodinger each E n (X) is a single- valued function of A but the critical points make 
the Richardson ones multi- valued functions of E. The Schrodinger critical points with 
their horizontal tangents have become square roots with their vertical tangents and 
singular derivatives. Since the oscillation counts are invariant under spectral inversion, 
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Figure 4. Positions of critical points for A > and n = 1 — 6. The two critical 
points that lie on the same real locus are connected by arrows. 

these counts in effect rotate with the graph and remain unchanged. 

In the Schrodinger problem it was argued that there were complex paths that 
allowed passage from one eigencurve to another and we expect the same type of 
connections here. As E becomes large and negative there will be two couplings ±A 
for the same E. This double-valuedness can be resolved by putting each branch on a 
different sheet so that we extend our notation, \ n (E) — > X^(E), where n is assigned 
by virtue of the oscillation count (n — 1) when E < and ± indicates the sign of 
the coupling in that region. As E becomes positive, the inverted Schrodinger critical 
points become square roots that couple one sheet to another. In figure[5]it is seen that 
the functions X+(E) and X~(E) have a common branch point at E = n 2 (w 2 /4) where 
the two couplings vanish. These are the images of the Schrodinger problem at A = 0. 
As in the Schrodinger case we expect an infinite number of sheets. By counting sheets 
or singularities it seems likely that each of the Riemann surfaces has infinite genus, 
and speaking crudely, spectral inversion takes one of them and turns it inside out to 
form the other. 

To guide our results when the eigencouplings leave the real domain, we will 
use a theorem (5.8.2) of Zettl 6J that gives a bound on the number of complex 
eigencouplings. That the Richardson problem is right-indefinite is not sufficient to 
guarantee complex results. It is necessary that the operator implicit on the left of the 
equation f4} be negative-definite; such a situation is termed left-indefinite and there 
is a simple test for it. Call Richardson with weight w(x) = — sgn(a;) the w(x) problem 
and, if a related problem with weight has m negative eigencouplings, then the 

w(x) problem can have at most 2m non-real eigencouplings. The |w(a;)| problem is 
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Figure 5. \ n (E) versus E for n = 1 — 4. 

right-definite and reads 

-<tt' n {x) + [\ n {E) + E]<t> n {x) = Q (21) 
4>„(±1) = 0, 

where X n (E) is the eigencoupling of the problem and the solution is X n (E) = 

n 2 (7r 2 /4) — E. The first value of E that allows complex results is one just to the right 
of 7r 2 /4 which is as expected and is where the 1 + and 1~ curves meet at their branch 
point. Note that the theorem confirms that nothing interesting happens for E < 7r 2 /4. 
In the general case if E is just to the right of n 2 (7r 2 /4), the bound is 2n. 

We now show results in figure [5] for \f(E) which was originally studied in [?]. 
Rather than trying to plot complex functions, we only show whether the result is real 
(R), imaginary (iR), or complex (C). It is seen that there are short segments in E 
in which the three spectral types occur over and over and this continues past that 
shown in the graph. The boundaries between segments are square root singularities 
that arise from inversion of Schrodinger critical points. The integer below each real 
segment is the oscillation count in that interval. That this changes from one section 
to another was considered a puzzle in earlier work 0]; whenever the eigencoupling 
returns from its complex plane, the eigenfunction comes back with another node. 
There are two ways of clarifying the puzzle. The Richardson way is to recall that the 
eigencurves of figure [5] inherited the oscillation counts from the Schrodinger curves. 
The behavior of Xi(E) represents real results that would be visible on eigencurves 
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separated by non-real segments that apparently induce a change from one eigencurve 
to the next one on the right. If viewed only in the real domain, this would appear as 
jumping from one curve to the next, with a unit increase in oscillation count. In the 
Schrodinger view there is a more natural mechanism for this curve jumping and that 
is the process of level crossing. The Schrodinger image of the Richardson segments 
would appear on figure Q] and involve movement on the A's, -B's and C's so that level 
crossing must occur between the real segments and this causes the needed change in 
oscillation count. The Richardson jumping occurs horizontally in figure[5]with E real 
and increasing, while Schrodinger jumping occurs vertically in figure[T]with E n (X) real 
and increasing. In either case a jump of curves involves a unit increase in oscillation 
count. 

iR 
I R 

ll Y ▼* C R C R C R 

-I H Slj-e & H & j & H Q 

TT t t i 

1 2 3 4 

i 1 1 1 1 1 1 1 £- 

20 40 60 80 

Re E 



Figure 6. \J{E) versus E. The symbols denote real (R), imaginary(iR), and 
complex (C). The integers give the oscillation count. 

The eigencoupling \^(E) need not be shown on a separate graph since the plot 
would not differ from that seen in figure [H In the real region of E < 7r 2 /4, we have 
by construction Xi (E) = — Xi(E) and we find this to be true also in the complex 
domain and also for other n that have been studied. This is the realization of the idea 
discussed in section [5] that for A on A, B, or C, there will be two couplings of opposite 
sign that occur on two different Richardson sheets. 

We combine in figure [7] results for X^(E), X^E), and X^(E) together with the 
Zettl bound (ZB) on the number of complex results. To the left of 7r 2 /4, the bound 
is zero and all results are real. In addition to X+(E) for n = 1 — 3, there is a X~(E) 
that is not shown and has opposite sign. In the first imaginary (iR) segment the two 
couplings Xf (E) saturate the bound at 2 and all other results are real. At E just above 
7r 2 , everything turns real and the bound is not achieved for a short interval. At E just 
above 57r 2 /4, the 1^ and 2^ segments become complex with the 1^ contributing ±A 
and the 2 ± giving ±A* and the bound of 4 is reached. Since the Schrodinger image of 
these intervals all involve n = 2 and 3, they are the Richardson version of the points 
singled out on the real loci of type C of figure [3] where we anticipated four couplings 
for some real E. This scheme continues on as E increases and more couplings become 
complex. 

It is useful to compare the complexity of the Richardson spectrum with what 
would be observed in the Schrodinger problem. There as real A varies we will have 
one spectral type (R), no change in oscillation count, and no square roots. The 
Schrodinger problem appears to be much simpler but our overall task was to show 
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Figure 7. The spectral types for Xn(E) for n = 1 — 3 versus B. ZB is the Zettl 
bound. 



that the two functions are inverses. This seeming paradox is avoided by observing 
that the quoted Schrodinger simplicity is all in the domain of the real loci of type A. 
It will turn out that Schrodinger is sufficiently complicated on the B's and C"s so that 
no conflict exists. 



7. Spectral inversion of Richardson 

We have argued above that the Schrodinger A's, B's, and C's with their real values 
and complicated places become the Richardson spectrum after inversion with its real 
places and complicated values. This argument does not specify which real locus is the 
image of a certain Richardson segment. The details of this will be worked out below 
in a few cases. 

Suppose we could go to the Schrodinger problem and erase the B's and C's from 
existence so that only the A's remain. The modified Richardson problem would then 
have only real eigencouplings and the system would lose much of its interest. We 
conclude that the following two statements arc equivalent: the Schrodinger problem 
has B's and C's and the spectrum of Richardson is complicated. This is consistent 
with our general argument that the two types of complication are the inversions of 
each other. 

The first four segments of the Richardson spectrum in figure [6] are K, iM., K, C 
with the oscillation count being and then 1 in the real segments and with the segment 
boundaries being square roots. In figure [5] are shown the n = 1 and n = 2 Riemann 
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sheets on which E n {\) takes its values. These sheets have a pair of common branch 
points on the imaginary axis while sheet n = 2 has four branch points leading to sheet 
n = 3. Every double-x locates a critical point. To match the Richardson sequence of 
values we need a Schrodinger sequence of places A G A, B, A. C, with the boundaries 
being critical points, the oscillation count changing from to 1, and with E n {\) always 
increasing. The point labeled E = on the positive real axis of the n = 1 sheet is 
the image of Xt{E) at E = 0. If we follow the arrows on the two sheets, we find the 
sequence of places A G A, B, A, C, the boundaries are critical points, and E n (X) is 
always increasing. When the branch point is reached on ImA < on the n = 1 sheet, 
the first instance of level crossing occurs where we do not encircle a branch point, we 
simply touch it, change sheets, E n (X) continues to increase and when the origin is 
reached on sheet n = 2, the oscillation count will be unity as required. The origin is 
a minimum so that we can go left until a maximum is reached. This maximum is a 
critical point so that a real locus of type C is available and the image point moves 
onto it and £2 (A) remains real and increasing. This completes the required sequence 
\ <E A, B, A, C. 



(a) n=l sheet 
Im X 




X 
V 



( h ) n = 2 sheet 
Im X 



/ 



"^K Re X^ 



I 
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Figure 8. Riemann surface of E n (X). 
broad lines are branch cuts. 



(a) n = 1 sheet, (b) n = 2 sheet. The 



Going back to the origin of the n = 2 sheet, the Schrodinger sequence A G 
A, C, A, C, ■ ■ ■ must be achieved to match the Richardson result of figure[Sl The way 
that this is accomplished is illustrated schematically in figure[9](a) for the n = 2 sheet 
and in figure [9] (b) for n = 3. The origin for n — 2 is a minimum, we proceed left to 
the maximum, move onto the real locus of type C, the path encloses a branch point of 
type n — 2 — 3, we pass to sheet n = 3, and when we return to the real axis in figure [9] 
(b), the oscillation count will be 2 and we will again be at a minimum. This sequence 
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of steps is repeated over and over leading to larger n and oscillation count. This is our 
qualitative argument that the two functions are inverses of each other. Other X^(E) 
could be studied in a similar way. 

(a) n=2 sheet (b) n=3 sheet 

Iffl X Im X 




Re X 



Figure 9. Schematic Schrodinger image of the Richardson function X^(E). (a) 
n = 2 sheet, (b) n = 3 sheet. 

The Schrodinger sequence of the places A S A, B, A, C, A, C, • ■ • is the 
realization of the curve jumping discussed in section [6l When viewed only in the 
real domain in figure [2 we would see the sequence A(n — 1), A(n — 2), A(n = 3) 
with the invisible level crossing sequence B, C, C providing the mechanism for getting 
from one eigencurve to the next. In the Richardson problem there is nothing as specific 
as level crossing to account for the change in oscillation count so that we rely on the 
Schrodinger description and argue that the count is invariant under spectral inversion. 

Here we add few additional comments. A real Richardson segment may not 
contain a critical point since that would imply that the Schrodinger problem is singular 
for real A. On an imaginary Richardson segment there is no such restriction and 
Xf (E) must have a critical point that is the image of the branch point in figure [8] 
(a) for ImA < 0. This critical point differs from the Schrodinger examples since two 
imaginary loci on which Kc\^(E) = 0, intersect orthogonally. In this way the local 
information on the two Schrodinger sheets near the square root gets unwound and 
inverted and appears on one Richardson sheet. On the neighboring sheet containing 
Ai (E) = -Xf(E), there is a similar critical point at the same value of E that is the 
image of the branch point for ImA > 0. The Schrodinger image of the function \i{E) 
starts on the left at E = on the n = 1 sheet and always moves in directions opposite 
to those of the arrows. 

It is amusing to conjecture that Xf (E) continues its oscillation R, C, R, C, • all 
the way to E = oo. The Schrodinger image of this would be an oscillation A, C, A, 
C, • all the way to n = oo and Re A = — oo. In this way the information on the real axis 
of one Richardson sheet would appear on an infinite number of Schrodinger sheets. If 
we single out the point where £a(A) = on the right side of the n = 1 sheet and also 
the point A = — oo on the asymptotic sheet, then these points would be connected by 
a very long real locus that is the image of positive Richardson real line. 
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8. Oscillation counts and defective multiplicity 

This section focuses on the observation - quite familiar by now - that the Richardson 
problem is more complicated for real E than the Schrodinger one is for real A. To 
examine this further, two questions are posed. 

(i) For given real E, how many Richardson eigenfunctions have oscillation count to? 

(ii) How many coincident square root singularities of the eigencoupling are there over 
a given real value of El 

If we ask these questions of Schrodinger quantities for a given real A, the very simple 
answers are: one for all to, and zero. Figures [1] and O contrast the singled-valued 
simplicity of Schrodinger with the multi-valued complexity of Richardson. 

We have no general result for the first question but we will work through some 
particular cases. For some real E let N m be the number of Richardson eigenfunctions 
with oscillation count to. In figure [5] if we approach a vertical tangent from the side 
with two real branches, there will be two eigenfunctions with the same count which 
limits to a single one at the branch point since the geometric multiplicity is unity. 
Some special cases are given in the following list that can be followed in figure [5] 



E < tt 2 /4 ps 2.47 




N m = 


-- 2, to > 




E = 7T 2 /4 




No = 


l;JV m = 2,m> 1 




tt 2 /4 < E < ir 2 ps 


9.87 


No = 


0; N n = 2, m > 1 




E = TT 2 




N = 


0;iVi = 3;JV m = 2,m 


> 2 


7T 2 < E < 5tt 2 /4 f 


« 12.34 


N = 


0;JVi = 4; N m = 2, to 


> 2 


E = 57r 2 /4 




N = 


0;JVi = 2;N m =2,m 


> 2 



We stop here because further cases involve numerical estimates for minima. Notice 
that N\ goes through the sequence 2,3,4,2 and iVi = after this. Things get very 
complicated as E increases and obtaining a general result seems difficult. 

In the Schrodinger problem it is not a significant event if there are two distinct 
couplings that are critical points with the same energy. Many cases can be seen 
in figure [T] and they arise from E n (— A) = E n (\) at critical points. Such an event 
when viewed in the Richardson problem is more significant mathematically since the 
eigencoupling will have coincident branch points on different sheets. The task posed 
by the second question is to count this number for given E. In dealing with this it is 
more convenient to use the Schrodinger description at critical points rather than the 
Richardson one at square roots. In problems of this sort these critical point positions 
and values are usually undecipherable numbers that are presumed irrational and can 
only be estimated numerically. The existence of the Binding and Volkmer [3] lattice of 
maxima means that we can reduce this part of the problem to integers, after dividing 
out a common irrationality. Apart from those at A = 0, we have no model for the 
minima and their positions and values are found to be in the undecipherable class. It 
follows that only about half of the problem can be treated exactly. Some critical point 
positions and values for minima with A 7^ are listed in table [5J 

In dealing with the second question we call the number of coincident square 
roots the defective multiplicity N4. Note that over almost all E, Nd — because no 
square roots are present. The first unexpected case that might be called "accidental" 
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Table 2. Position of minima of E n (X) for A ^ 0. 



n 


A 


E n {\) 


3 


±6.151546 


21.996039 


4 


±23.271272 


41.909383 


5 


±10.258305 


61.478860 


5 


±52.084097 


71.536568 


6 


±34.746014 


91.264186 


6 


±91.259508 


111.019851 



occurs when the energy at the Schrodinger outer maxima for n = 4 at A = ±6tt 2 
coincides with that for the central maximum of -Es(Q) = 25(7r 2 /4), leading to Nd = 3 
in Richardson. Because of symmetry, non-vanishing Nd will be even and at least 2 
unless a value at A = is involved, as it is here. 

To study further examples we set up the lattice of maxima shown as dots in figure 
[Tniwith the minima shown as crosses. Lattice points are located by the pair of positive 
integers (i,j) giving the couplings A^-, the energies Eij and the quantum numbers mj 
as follows 

\ ij = (ir 2 /4)[2(i + j-l)(j-i)}, 

= {tt 2 /A)[2i 2 + 2j 2 -2i- 2j + 1), (22) 
mj = i + j — l. 

Dividing out the irrationality gives integer values for the reduced couplings and ener- 
gies 

Xij = 2(i + j — — i) = even integer, (23) 
Eij = 2i 2 + 2j 2 - 2i - 2j + 1 = odd integer. (24) 



The defective multiplicity is the number of distinct pairs of positive integers (i,j) that 
have the same E^. We can form a symmetric matrix E whose first few elements are 



E 



1 5 13 25 

5 9 17 29 

13 17 25 37 

25 29 37 49 



(25) 



Note that Nd — for any real number that is not present in E, at least as far 
as maxima are concerned. By inspection, Nd = 1 for Eij = 1,9,49; Nd = 2 for 
Eij = 5, 13, 17, 29, 37; and Nd = 3 for Eij = 25, which is the case we called accidental. 
In general we must have Nd > 2 for i ^ j and Nd > 1 and odd if i = j. 

A classical problem in number theory asks how many positive, negative or 
vanishing integer pairs lead to the same N = i 2 + j 2 . The well-known solution 
[14] involves determining all the integer factors that divide N with no remainder. The 
even divisors are ignored and the odd ones are separated into two classes 
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Figure 10. Positions of extrema of E n (X) versus A. The solid dots are maxima 
and the crosses are minima. The lattice coordinates arc indexed as 

rn(N) = 1 mod 4 = 1,5,9,13, 

n 3 (N) =3 mod 4 = 3, 7, 11, 15.... (26) 

The classical result for the number of pairs is 4[ni(N) — 71,3 (TV)]. Our problem dif- 
fers because only positive integers are involved and our expression for Eij is a circle 
centered at — (1/2,1/2) rather than one centered at the origin. By studying 

hundreds of examples we find empirically that a similar result is obtained 

N d = n 1 (E ij )-n 3 (E ij ). (27) 

For example, the factors of E^ — 25 are 1, 5 and 25 which are all congruent to 1 mod 4 
so that Nd = 3, as has already been seen. If we return to the matrix E and choose 
an odd integer that is not present such as 15, we obtain m(15) = 713(15) = 2 so that 
N d — is verified since the energy is not on the lattice. 

As we consider larger E^ , the first appearance of a certain Nd is found to grow 
slowly. Table [3] gives the smallest value of Eij that produces the listed values of 
Nd and the first three entries have already been discussed. It appears that N d goes 
slowly to infinity with Eij. The last entry with Nd = 10 involves maxima with 
coordinates (i, j) = (6, 64), (13, 63), (28, 58), (30, 37), (43, 48), and the number doubles 
by exchanging factors. The quantum numbers involved are n = 69, 75, 85, 86 and 90. 
The total number of maxima for the eigencurves is 405 and of these, 10 are coincident 
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Table 3. Smallest Eij that has defective multiplicity Nd- 



Eij 


N d 




N d 


1 


1 


325 


6 


5 


2 


15625 


7 


25 


3 


1105 


8 


65 


4 


4225 


9 


625 


5 


8125 


10 



at = 8125. These curves must oscillate in a very complicated mariner so that, 
although they enter each others vicinity, they must do so without intersection. 

To this point we have only studied N d that arise from coincident Schrodinger 
maxima of the form (max, max) but we have ignored the types (max, min) or (min, 
min). The only exact information we have about minima is that at A = 0, the reduced 
energies are even integers with £^(0) = n 2 , with n = 2,4, 6 • • ■. Since at maxima 
E.^ = odd integer, (min, max) degeneracy of this type is not possible. For A ^ 0, a 
(min, max) degeneracy requires E(X) — odd integer at the minimum. For n < 10, we 
find that this does not happen and results at the minima are all of the undecipherable 
type. There will be many degeneracies of the (min, min) type due to E n (—X) = E n (A). 
To find others would involve problems with numerical round-off error. We are confident 
that there are none of this type for n < 10. 

9. Summary 

In this paper we have studied the eigenvalues of two SL problems, each dependent 
on a parameter, and shown that they are inverse functions of each other. This has 
allowed an understanding of the mysterious features of the Richardson spectrum in 
terms of Schrodinger behavior in the complex domain of the coupling parameter. In 
Schrodinger the properties of E n (X) on the three types of real locus were shown to 
be related to the three spectral types of Richardson. The mathematical features that 
generate the complexities in each description are the critical points and square root 
branch points of the eigenquantites. In Schrodinger the critical points act as sources of 
the B's and C's while the branch points tie the Riemann sheets together and facilitate 
the level crossing that accounts for the oscillation count variability of the Richardson 
spectrum. The spectral inversion of the Schrodinger critical points yields Richardson 
square roots that separate the three spectral types and tie those Riemann sheets 
together 

The complications in each problem are driven by the simple feature that the 
Richardson weight, or equivalently the Schrodinger potential, changes sign once. We 
saw in section @] that these complexities were absent for problems with no sign change 

If asked to explain the Richardson spectrum on its own terms and with no 
reference to another problem, we do not know how to proceed. It may be best to 
think of the two problems as a closely related pair, with Schrodinger contributing to 
the understanding of Richardson. 
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